Slow-light effect in dual-periodic photonic lattice 
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We present analytical and numerical study of photonic lattice with short- and long-range 
harmonic modulations of the refractive index. Such structures can be prepared experimentally 
with holographic photolithography. In the spectral region of photonic bandgap of the underlying 
single-periodic crystal, we observe a series of bands with anomalously small dispersion. The related 
slow-light effect is attributed to the long-range modulation in the photonic lattice, that leads to 
formation of an array of evanescently-coupled high-Q cavities. The band structure of the lattice is 
studied with several techniques: (i) transfer matrix approach; (ii) analysis of resonant coupling in 
process of band folding; (iii) effective medium approach based on coupled-mode theory; and (iv) 
Bogolyubov-Mitropolsky approach. The latter method, commonly used in the studies of nonlinear 
oscillators, was employed to investigate the behavior of the eigen-function envelopes and the band 
structure of our dual-periodic photonic lattice. We show that reliable results can be obtained even 
in the case of large refractive index modulation. 
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INTRODUCTION 

One of the attractive features of photonic crystals 
(PhCs) is the possibility to alter the dispersion of the 
electromagnetic waves [Ij. "Slow light" [3] is the property 
that attracted a great deal of practical interest as it leads 
to low-threshold lasingjj, |j, Q, pulse delay [6|, |7|i opti- 
cal memories [8|, and to enhanced nonlinear interactions 
[Sl, llO| . There are several approaches to slowing down the 
light in PhC structures: 

(i) At the edge of a bandgap, the dispersion (group 
velocity) Vg = dco/dK becomes zero. This property has 
been extensively studied and used in practice for control 
of spontaneous emission llj and gain enhancement in 
lasers [3J, y, |5| . In the vicinity of the band-edge the group 
velocity strongly depends on frequency leading to strong 
distortion of a pulsed signal, rendering it impractical for 
e.g. information processing applications. 

(ii) Unlike one-dimensional PhC, in two- and three- 
dimensional structures the dispersion of the high- 
frequency bands becomes smaller over the finite region 
of frequencies. This is due to the increased density of 
states which leads to band-interaction and flattening of 
the photonic bands [1^, ll3| . These high-frequency fea- 
tures, however, are of limited use in practice because they 
allow little control, are not spectrally isolated, and are 
usually susceptible to disorder [l4|. 

(iii) Based on the Coupled Resonator Optical Waveg- 
uide (CROW)[i El, 111113 idea, the periodic array of 



weakly coupled optical resonators made in a photonic 
crvstal[18lll9l.l2(J| mitigates the limitations in approaches 
(i-ii). Hybridization of high-Q resonances due to periodi- 
cally positioned structural defects can form a defect-band 
in the spectral region of the photonic bandgap: 



uj{q) = rj[l + KCOs{qa)] 



(1) 



As a result, the group velocity is reduced over a substan- 
tial range of frequencies, whereas the band characteristics 
can be controlled via the parameters of the structural de- 
fects. 

A periodic arrangement of structural defects in the 
photonic crystal creates a photonic super-crystal (PhSC) 
with short-range quasi-periodicity on the scale of the lat- 
tice constant and long-ran ge p eriodicity on the defect 



separation scale [2l|, 122, l23|, l2J, l25|, iM l27| . We recently 
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proposed PhSC with dual-harmonic modulation of re- 
fractive index[23| that can be fabricated[29| via a single- 
step interference photolithography technique. Because 
all resonators are produced at once, our design minimizes 
fabrication error margin and ensures the large-scale pe- 
riodicity which is essential for hybridization of the reso- 
nances of individual cavities in an experiment. 

The purpose of this work is to theoretically investigate 
the optical properties of one-dimensional PhSC with a 
combination of analytical and numerical techniques. We 
consider systems where the dielectric function is given by 

e(x) = £o ^ — [1 + 7C0S (2ttx/L)] [1 + cos (2TTx/a)] . 

1+7 

(2) 
Here Eq is the background dielectric constant. The ampli- 
tude of the short-range (on scale a) modulation gradually 
changes from Ae x (1 — 7)/(l + 7) to Ae. L sets the scale 
of the long-range modulation and we will concentrate on 



L/a ^ 1 limit. An interesting situation may arise when 
L/a is non-integer. The spectrum of such structures ex- 
hibits the characteristic "butterfly" structure [23J, similar 
to the spectrum of an electron in a solid with an applied 
magnetic field [30|. We, however, will limit our consider- 
ation to commensurate lattices where L/a = iV is a large 
integer. 

The paper is organized as follows. In Section 2 we dis- 
cuss the possibility of reproducing the proposed struc- 
tures in holographic photolithography experiments. In 
Section 3 we investigate the optical spectra of finite and 
the photonic band structure of the infinite PhSC with 
transfer matrix technique. In Section 4 we study pho- 
tonic band folding for moderate values of N and discuss 
the onset of the formation of photonic bands with small 
dispersion. In Section 5 we apply perturbative, "effec- 
tive medium" , approximation to our dual periodic crys- 
tals with large N . Using separation of scales L ^ a 
technique from nonlinear oscillator theory, in Section 6 
we derive a set of first order differential equations that 
describe our problem on the language of eigen-function 
envelope. We conclude with the summary Section 7. 



^1,2; (fi) long-range modulation L via 9i — 62', (iii) depth 
of the long-range modulation via relative intensity of the 
beams Ei/E2- Although the exact e{x) dependence in 
Eq. m differs from that in Eq. [21 it shows the same spec- 
tral composition and modulation property. If the above 
structure is realized in experiment, the discrepancy is ex- 
pected to cause only small deviation from the analytical 
results obtained further in this work, furthermore, it also 
becomes insignificant in the limit L ^ a. From now on, 
for the reasons of computational convenience we will only 
consider dielectric function in form of Eq. [51 

Dual-periodic harmonic modulation of refractive index 
can be also experimentally realized in optically-induced 
photorefractive crystals [3 Ij. Although, unhke QDPL[29|, 
the obtained index contrast is several orders of magnitude 
less, the superlattices created in the photorefractive ma- 
terials offer the possibility of dynamical control. Study 
of dynamical and nonlinear phenomena in dual-periodic 
lattices presents significant interest. However, it goes be- 
yond the scope of present publication and will not be 
considered in this work. 



AN EXPERIMENTAL REALIZATION OF PHSC 

In our previous work [28] we proposed to employ holo- 
graphic photolithography to fabricate a modulated pho- 
tonic crystal, PhSC. The experimental quantum dot UV 
photolithography (QDPL) technique - a bottom-up tech- 
nique that allows to synthesize quantum dots and other 
nanoparticles in selected regions of porous matrices [29J - 
is suitable for this purpose. The change of the dielectric 
function is related to the local density of quantum dots 
and it is proportional to the local intensity of the elec- 
tromagnetic wave used when exposing the sample. We 
studied a setup where four S-polarized laser beams are 
defined by 



<iL2,EL2 
Qfil , Eri 
qfl2 , Eii2 



ko{-sm{ei),0,cos{ei)},Ei 

A:o{-sin(02),O,cos(^2)},^2 
A:o{sin(0i),O,cos(0i)},Si 

A:o{sin(02),O,cos(02)},S2 



(3) 



Here q and E are the k-vector and amplitude of 
the beams respectively. Their interference Etot{x) ex: 
acos{kix) + P cos{k2x) leads to the following expression 
for dielectric constant 



e{x) = £0 + Ae [a cos(A:ia;) -f /3 cos(fc2a;)]^ 



(4) 



where a = Ei/{Ei + E2), /3 = E2/{Ei + E2), fci = 
ko sin(0i) and fc2 = /cq sin(02)- Choosing ki^2 — Tr/ain/L 
and a = 7/ [1 + 7] creates PhSC with dual-harmonic 
modulations, short- and long-range, similar to that in 
Eq. [21 Manipulation of the lithographic beams allows 
for easy control over the structural properties of the re- 
sultant PhC: (i) fundamental periodicity a via fcg and 



TRANSFER MATRIX ANALYSIS 

Transmission/reflection spectra of PhSC of finite 
length as well as the band structure of its infinite coun- 
terpart can be obtained numerically via transfer matrix 
approach. The propagation of the field through a in- 
finitesimal segment of length dx is described by [33 



M{x, x+dx) 



cos{kn(x)dx) n ^(x)sin{kn(x)dx) 
—n(x) sm{kn{x)dx) cos{kn{x)dx) 

(5) 
where we assumed that refractive index n(x) does not 
change appreciably over that distance. The matrix 
M{x,x + dx) relates the electric field and its spacial 
derivative {E, 1/k dE/dx} at x -\- dx and x. Total trans- 
fer matrix of the finite system is given by the product of 
the individual matrices 



Mt. 



= Y[ M{x, X + dx) 



(6) 



x=0 



Since in our case refractive index n{x) = e 



[IJa), is not a piece- wise constant (as in Ref. |25l. 
a continuous function of coordinate, one has to resort 
to numerical simulations. In what follows we apply ei- 
ther scattering or periodic boundary conditions to obtain 
transmission coefficient and Bloch number (7(0;) respec- 
tively. 

Fig. [2ta) plots transmission coefficient through one pe- 
riod of a dual-periodic system shown in Fig. [H A series 
of progressively sharper resonances occur on the lower 
or upper edge of the spectral gap of underlying single- 
periodic structure. Whether the peaks occur at lower or 
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FIG. 1: (a) Dependence of the index of refraction in dual- 
periodic photonic crystal as defined by Eq. (2). We used 
eo = 2.25, Ae = 1, A'' = 80 and the modulation parameter 
7 is equal to 0.25. (b) Local (position-dependent) photonic 
bandgap diagram for n(x) in (a). A- and B^ mark the 
frequencies of the foremost photonic bands on long- and short- 
wavelength sides of PBG of the corresponding single-periodic 
crystal. 
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FIG. 2: (a) Transmission coefEcSmt through a finite segment 
of length L (a period) of the periodic super-structure defined 
in Fig. [l] Solid and dashed lines corresponds to < a; < Na 
and —Na/2 < x < Na/2 segments (shown in the inset of 
panel (b)) respectively, (b) Solid and dashed thin lines plot 
the corresponding phase of t{u}). Bold line depicts the Bloch 
number g(a;) x a of the infinite crystal computed using Eq. [T] 



upper band edge depends on the particular definition of 
the unit cell as shown in the inset of Fig. [UJa). One 
can obtain an insight into this effect by examining the 
modulation of the spectral position of the "local" pho- 
tonic bandgap with coordinate as shown in Fig. [Ijb). 
This analysis is meaningful on a ^ Aa; <C L = Na scale. 
This condition can be satisfied in our case of slow mod- 
ulation - large N. At the frequencies such as Ai in 
Fig. m^b) the wave propagation is allowed in the vicin- 
ity of X — aN X (1/2 + m), whereas x — aN x m are 
locally forbidden, m being an integer. When considering 
< a; < Na segment of the lattice, resonant tunneling 
via electromagnetic states of the cavity at the geometri- 
cal center, A^ , leads to the low-frequency peaks in the 
transmission coefficient - solid line in Fig. [2lja).0n the 
other hand, transmission through —Na/2 < x < Na/2 
segment exhibits a series of sharp resonances which cor- 
respond to tunneling via B^ cavity states in the high- 
frequency region. 

The transmission coefficient through a finite segment 
of length L (equal to one period) can be related to the 
band structure of the corresponding periodic lattice |26| 



cos{q{uj)L) — Re\ 



1 



1 



tiio)' |t(^)| 



cos(0(a;)), (7) 



where we introduced the phase of the transmission coeffi- 
cient through t — \t\ exp[i(/<]. Fig. [3]shows that hybridiza- 
tion of the cavity resonances considered above leads to 
formation of the flat photonic bands. Their low disper- 
sion, small group velocity, may be exploited [28[ for prac- 
tical applications. 

In the vicinity of an isolated transmission resonance, 
t(w) is given as the Lorentzian 



tiu) 



i)^(r/2) 



i{T/2)-{lo-ujo) 



(8) 



where F is the full width at half maximum (FWHM) of 
the resonance and loq is the resonant frequency. Substi- 
tution of Eq. [8] into Eq. [7] gives the flat band described 
by 



tt'(g) = cJo[l ± KCOs{qaL)] 



(9) 



F 1 

where k — = 7^ "^ 1; ^-^d Q is the cavity Q-factor. 

Thus, the decrease of group velocity in the PhSC is di- 
rectly related to the increase of confinement and decrease 
of the coupling between the neighboring cavities. In our 
PhSC both these factors are described by the same pa- 
rameter - the cavity Q-factor. In a single-periodic PhC of 
finite length, the Q-factor of a band-edge mode depends 



on the system size. Comparing Figs. [lJb)l3]one can see 

that Ai ,Bl ' modes are in fact band edge modes in 
their intervals of free propagation. In our case L gives 
the characteristic length, which, as we will demonstrate 
below, also determines the mode frequency. As N in- 
creases, the eigenfrequencies of the modes shift towards 
bandgap. The associated decrease of the local group ve- 
locity contributes to the increase of the Q-factor of the 
resonators and further reduction of the group velocity in 



IE(x)r 



A, 



[N) 



B?^^ bands in N 



oo limit. 



Eq. [7] suggests that the dispersion relation Lij{q) is in- 
dependent of how segment of length L (the period of our 
structure) is chosen. However, the transmission coeffi- 
cient through < a; < A^o and -Na/2 < x < Na/2 
segments of the crystal show very different spectral com- 
position, Fig. Ufa). In order to understand how these 
markedly different fmictions lead the same uj{q), we will 
analyze the phase of the transmission coefficient 0, Fig. 

Hb). 

In a one-dimensional periodic system such as ours, the 
wave number qiuj) in Eq. [7] is equal to the integrated 
density of electromagnetic states. It is by definition a 
monotonically increasing function of the frequency in ex- 
tended Brillouin zone scheme. In PhSC, q x L increases 
by TT every time the frequency is increased through an 
allowed band, c.f. bold line in Fig. ^h). At the fre- 
quency in the middle of the band cos{qL) — because 
qx L = TT X {771+ 1/2). From Eq. ([7]) one can see, that (f> 
should be equal to tt x (rn + 1/2) at the same frequency. 
In the finite system mode counting phase defined 33] as 
tan((/)) — E' / E coincides with the phase of the trans- 
mission coefficient = that explains the monotonic 
behavior of 0(w). Eq. [7] leads to the fact that the quasi- 
states of the finite system occur at the same place as the 
corresponding band center of the lattice, irrespective of 
the definition of the unit cell. Thus, as it can be also seen 
from Fig. [SJb), 0(ijj) and q{uj)L intersect at vrx (m+1/2). 

Taylor expansion of the phase around frequency at the 
center of a pass band, wq, where q{uj)L = n x (m + 1/2) 
gives 



cos{q{uj)L) — (u; — LUo) X 



|i(ijJo)| duj 



(10) 



Here, the term that contained d|i(ti;o)| /dto dropped out 
because cos{q{uJo)L) — 0. Comparing Eqs. (|9ll0p shows 
that |i(cL;o)|~ d(f>{uJo) / duj determines Q = l/n and not 
only |t(a;o)|- Suppressed transmission compensates for 
slow phase change {e.g. solid line in Fig. [D(b) in the 
high frequency spectral region) and leads to the identical 
q{u!) for two different definitions of the unit cell. 

We also note that if the segment is chosen such that 
the correspondent "cavity" is located in the geometrical 
center {\t{uio)\ = 1), FWHM of the resonance (F) in the 
transmission coefficient is equal to the width of the pass 
band in the periodic lattice. This fact follows from Eqs. 
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FIG. 3: The left panel shows dispersion of PhSC ijj{q) reduced 
to the first Brillouin zone. The eigen-niodes which correspond 
to the series of flat bands in the vicinity of the parent-bandgap 
of the single periodic crystal are depicted on the right. Cal- 
culations were performed for the structure described in Fig. 

HI 



(|7l8p . It further emphasizes the analogy with CROW 
structures that we explored in Ref. ,2S]. 

We conclude this section by noting that long-range 
refractive index modulation creates alternating spatial 
regions that can serve as resonators separated by the 
tunneling barriers. Hybridization of the cavity reso- 
nances creates a series of photonic bands with low dis- 
persion. The envelope of the eigen-states in these bands 
Al \b\ ' is a slowly varying function of the coordinate, 
see Fig. [31 This effect stems from states proximity to 
photonic bandedge of the underlying single-periodic lat- 
tice. Possibility of separation of length scales into short 
(a of rapid field oscillations) and long [L of the slow am- 
plitude variation) will inform our further analytical stud- 
ies presented in the following sections. 



RESONANT APPROXIMATION 

Forbidden gaps in spectra of a periodic system arise 
due to resonant interaction of the wave with its Bragg- 
scattered counterpart [3J . The scattered wave appears 
due to presence of Fourier harmonic in the spectrum of 
the periodic "potential" , which is in the case of Helmholtz 



Resonant q, (even A'^) 
(odd N) 
Coupled components 
Coupling harmonics 
Center frequency, uq 



E- 



TV/L TV/L 

n/L 

,Em-{s+1) Es-N,Es -E_(s + i),i5jv- 



e_(jv_i), Eiv-i £-N, £n e_(jv+i), en+i 
\/eL \/eL \/eL 



Normalized width, Alu/ljo £(jv_i)/e £jv/e 



-(iV+l) 



/£ 



TABLE I: Results of resonant approximation analysis of Eq. [14] with dielectric function given by Eq. [2] Three columns 
correspond to the three resonant photonic band gaps that appear in the spectrum of the dual-periodic PhSC. The expressions 
hold for both even and odd A'' for the choice of parameter s: N — 2s and A'^ = 2s + 1 respectively. 



equation 



E"{x) + -^de{x)E{x) = -^sE{x), 



(11) 



is represented by {cj^ / 0^)66 {x) = {uP^jc^) \e{x^ — e]. Here 
we introduced the average value of the dielectric function 
e = e(x) — eo + Ae/[2(l-|-7)]. When 7 = 0, the condition 
Ae/e ^ 1 is sufficient to obtain the position and width 
of spectral gaps. Otherwise, an additional condition N x 
Ae/e ^ 1 needs to be satisfied instead. We will discuss 
the physical meaning of this condition at the end of this 
section. 

We begin by noticing that 6e(x) of our choice (Eq. ^ 
contains 



e{x) 



E 



£m exp 



.27r 

i——mx 
Ij 



(12) 



only eight nonzero Fourier harmonics: m — {±1, ±{N — 
1), ±N, ±(iV + 1))}. This fact allows one to do the ex- 



J 



haustive study of all resonant interactions as follows. Ex- 
pressing E{x) in terms of its Fourier components 



E[x) — eyi'p[iq{uj)x\ 2, ^m exp 



27r 
i — mx 
1^ 



(13) 



leads to an infinite system of linear coupled equations 



q{uj) 



27r 
— 1 
L 



En 



LO 






■)/ -Cv7To,_7- 



(14) 

where q is the Bloch number that varies in the first Bril- 
louin zone [0, tt/L]. For the extreme values of q there 
exists spectral range where the term in brackets of Eq. 
1141 can become simultaneously small for certain m and 
—m at q = and for m and —m — 1 at q = vr/L. If 
e{x) contains harmonic Sm' such that it couples these 
two Fourier components, the overall infinite system Eq. 
[14] can be reduced to two resonant equations. 



The results of such analysis are summarized in Table |T] 
and the correspondent band structure is shown in Fig. [4] 
Introduction of the long range modulation in the dielec- 
tric constant results in expansion of the unit cell from a 
to L — Na and, thus, the reduction of the Brillouin zone 
accompanied by the folding of photonic bands. The cases 
of even N — 2s and odd A^ = 2s -|- 1 should be distin- 
guished. In the former, the primary photonic bandgap 
(lie) of the single-periodic lattice is reappears a,t q = 0, 
whereas in the latter (IIq) it is located at q — tt/L. Our 
analysis shows that the nearest in frequency gaps, Io,o 
and IIIo,o, also become resonant. Fig. [3] For the consid- 
ered refractive index modulation of Eq. [51 the normalized 

7Ae/4£ 



width of the satellite gaps is smaller e(jv-i)/£ ■ 



1 + 7 



than that of the central gap by a factor of 7. This pa- 
rameter by definition is less than unity. 

We can see that folding and the offset of the forma- 



tion of flat A\^ t Bi ' bands is captured in this approx- 
imation. The criterion of its applicability can be found 
by considering the contributions of non resonant terms 
in Eq. [14] We find that for all three gaps the criteria 
are qualitatively the same. Therefore, we present the de- 
tailed analysis of only one particular resonance, IIIc. The 
condition that the closest non-resonant Fourier compo- 
nents E-s-2, E-s, Es-1 and E^+i are smaller than the 
resonant ones E^s-i and Eg leads to 



ANe- 



<1. 



(15) 



2{N + lfeN-i 
In the limit of very large N the second term in the de- 
nominator becomes dominant and this condition cannot 
be satisfied regardless of the value of Ae. Thus, N should 
be finite. It can be seen that the condition under which 
the first term in the denominator is dominant, is consis- 
tent with the entire inequality Eq. [15] and can be put to 



EFFECTIVE MEDIUM APPROXIMATION 
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FIG. 4: Dispersion relation computed with transfer matrix 
formalism for eo = 2.25, Ae = 0.32, N = 9 (dashed line) and 
N = 10 (solid line) and the modulation parameter 7 is equal 
to 0.25. For this set of parameters, the applicability condition 
Eq. [15] of the resonant approximation is satisfied. 



read Nsi/e ^ 1. Taking the case most restrictive £i we 
finally obtain 



Ae 



x7V<l, 



(16) 



where we neglected by 7 for simplicity. 



The derived expression Eq. [16] has a clear physical 
meaning. Indeed, from Table HI the frequency of the 
bandgaps I and III approach the central gap inherited 
from the single periodic system as 1/A^. At some point a 
bandgap of width Alo/uq = Ei/e begins to substantially 
perturb the pass band of width qmax x c ~ loq/N sepa- 
rating the consecutive gaps. The resonant approximation 
breaks when these two scales become comparable. This 
condition results in Eq. 1161 In other words, the approx- 
imation considered in this section can at most capture 
the onset of flattening trend in A]^ \ B\ bands and 
fails when N is increased to the degree when these states 
become abnormally flat - Auj/qmax ^ c/e throughout 
the band. More sophisticated approaches are considered 
below. 



Gratings written in the core of photosensitive optical 
fibers are often analyzed with the help of coupled-mode 
theory (CMT):35|. In both shallow gratings with long- 
range modulation in fibers[25!, i27|| and our PhSC, the 
forward and backward propagating (locally) waves con- 
tinuously scatter into each other. The advantage of GMT 
is in that it considers the amplitudes of the forward and 
backward waves directly. This tremendously simplifies 
Maxwell equations. Ref. [36[ also considered the fiber 
gratings with deep piece- wise constant index modulation. 
In this section we employ the GMT-based method devel- 
oped by Sipe et al [25| to obtain the spectral positions of 
the flat photonic bands formed in PhSG. 

For shallow modulation^ small Ae, our Eq. ^ can 
be br oug ht to resemble the model function considered in 
Ref. [ill 

n{x)/nQ = 1 + (t{x) + 2k{x) cos [2kox + ip{x)] (17) 

with the following choice of parameters 

7Ae/4 



cr(x) = 


1+7 2tt ^ 


As/2 " '"' L ^' 

1 + 7 




Ae/8 


k{x) = 


1+7 f. , 2^ 
x(^l+7Cos X 

1 + 7 


(fix) E 


/ Ac/2\'/' 
I 0; no - eo + _, , ; 



(18) 



fco — 7r/a. 



GMT of Ref. [25!| is applicable as long as these functions 
have slow dependence on x, on the scale much larger then 
k^^. This condition is indeed satisfied in PhSG with 

By introducing small detuning parameter 



UJ — UJQ 



<1, tJo 



kpc 
no 



we can, following Ref. [2^], obtain the governing equa- 
tion for the quantity -Ee// related to the envelope of the 
electric field 



^e// 



dx^ 



+ klnlff{x,i^)Eeff =0. 



(19) 



Frequency and position dependent effective refractive in- 
dex 



^e// 



= {{a{x) + Af - n{xfY'^ 



(20) 



determines whether propagation is locally allowed (real 
neff) or forbidden (imaginary n^ff). This is similar to 
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FIG. 5: (a) The value of integral in Eq. (|2ip . solid line, as 
a function of frequency is shown. For an easy comparison 
with (b), the plot is transposed so that w is plotted along 
y-axis. The open circles depict frequencies that satisfy the 
quantization condition of Eq. (|2ip . The dashed lines denote 
the actual position of the photonic states, as determined by 
direct numerical analysis of Section 3. (b) Color-coded plot 
of Relneff{x,u})] given by Eq (f20)) . The solid line shows the 
boundary of the region where Im[neff{x,uj)] =^ 0. For com- 
parison we also show the local PBG of Fig. [Ub), dashed line. 
In both (a) and (b), the parameters of Fig. [T]are adopted. 



ically. The open circles denote the the frequencies at 
which the quantization condition of Eq. (pT|) is satisfied. 
In the system with the parameters we used for illustration 
in previous sections, the obtained solutions are in good 
agreement with numerical results obtained with transfer 
matrix approach, Section 3. This suggests that index 
variation given by Ae = 1, Eq = 2.25 was sufficiently 
small for this approach to be applicable. 

We finish the current section by noting that it would 
be desirable to retain the attractive property of CMT en- 
velope approach without being constrained by the con- 
dition of small refractive index modulation. The latter 
may not always be justified in the experimental situation 
of interest [29|. In the following section we develop such 
an approach. 



BOGOLYUBOV-MITROPOLSKY APPROACH 

In this section we will consider the standing-wave so- 
lutions of Eq. (fTT|) . In this case, the corresponding 
E{x) can be chosen a real function by appropriate choice 
of normalization. Than, we make Bogolyubov anzatz 
MM: 



E{x) ~ A{x) cos {kox + (j){x)) 
dE{x)/dx — —kQA{x)sin{kQX + (l){x)). 



(22) 



our definition of local PBG diagram which we studied nu- 
merically in Section 3 Fig. [Sb compares CMT's region of 
evanescent propagation (solid lines) to the numerical cal- 
culation (dashed lines). We attribute the relatively small 
discrepancy to the assumption of shallow modulation we 
made in arriving to Eq. ([20)1 . 

Eq. (fT9|) is formally similar to Schrodinger equation. 
Our previous analysis shows that the single-period states 
associated with photonic bands A- \Bi ' are confined 
to the region of classically allowed propagation in the 
language of quantum mechanics. By analogy, WKB ap- 
proach of quantum mechanics can be applied ^2^] to 
determining quantization of energies inside our optical 
equivalent of a quantum well 



fco 



rieff (x, ijj)dx — {ra + l/2)7r 



(21) 



Xl and xr are respectively the left and right turning 
points defined by the condition neff{xL,R,^) = 0, to is 
an integer. The solid line in Fig. [5^ depicts the value 
of the integral, x-axis, as function of lo obtained numer- 



where as in the preceding section, Uq = ii/a. The above 
equations defines the amplitude and phase functions. Its 
substitution into Eq. ([T4|) gives so-called Bogolyubov 



equations in the standard form|37l. |^ 



d(j){x) 

dx 
dA{x) 

dx 



1 

ko 

2fcn 






A-2 



—e{x) 



cos^ (fcox -I- ^{x)) (23) 
sin 2 (k^x ~\- (t>{x)) , 



No approximations have been made so far. The structure 
of the above equation suggests that conditions dA/dx ^ 
koA and d(j)/dx <C kQ(f) can be satisfied in the vicinity of 

the spectral region where (l/ko) uj'^/c'^e{x) — k^ <^ fco- 
Overbar denotes the average over one period. Compari- 
son with the analysis in the previous sections shows that 
this condition is satisfied in the vicinity of the primary 
photonic bandgap. In the system of interest A^ :^ 1, 
this observation justifies "averaging-out" the fast spectral 
components - Mitropolsky technique [38|. The averag- 
ing procedure leads to the following system of nonlinear 
equations on the slow-varying amplitude and phase 



d4>{x) 
dx 
d\ogA{x) 
dx 



1 

2ko 



^-£0 - ^0 



1 w2Ae/2 



u? Ae/2 / 27r 

-^ 1 + 7 cos --a; 

27r \ . ,, , 
7 cos —a: I sin 2q)\x). 



- cos2(f){x) 



(24) 
(25) 



r 




FIG. 6: Filled circles in panels a,e denote the spacial position 
where the particular (j>{x) turns to rm:/2. At these special 
points dA{x)/dx = denoted by the vertical dashed lines in 
b-d and f-g panels. 



In derivation of Eqs. (|24l25p we used the explicit form of 
e(a;) given by Eq. ^. 

We begin the analysis of the obtained Eqs. (I24I25|) 
with a discussion of the appropriate boundary conditions. 
In deriving these equations we limited the consideration 
to the real- valued solutions of the original Eq. (fT4| , that 
can be found only for a discrete set of frequencies. At the 
sought frequencies, the corresponding amplitude function 
should incur the periodicity of the dielectric function Eq. 
^ that implies 



0(i) = 0(0) + rmr 
sin2<?!)(0,L/2,L) = 0. 



(26) 
(27) 



The first condition is obtained by requiring sin 2(j){x) in 
Eq. (|25p to be periodic. Symmetry of the modulation 
profile A(x), see Fig. [TJa), and continuity of its derivative 



lead to dA{x = 0, L/2, L)/dx = condition. This can 
only be satisfied by requiring Eq. (P7|) . because other 
factors on the right hand side of Eq. (pS)) are positive 
functions. 

The equation which determines evolution of the phase, 
Eq. (PH) . is self-contained. Hence, its solution together 
with the constraints given by Eqs. (|26l27p is sufficient to 
obtain the spectrum of system and (/'(a;). The amplitude 
is to be recovered in the second step by simple integration 
of Eq. ((25|l with the found phase 4>{x). 

Fig. E] shows the solutions of the Eqs. (|24l25|26|27p ob- 
tained by fourth order Runge-Kutta numerical method. 
In accord with our expectation, for each band there ex- 
ist two solutions 4'{x) which correspond to the standing- 
wave bandedge modes at fc = 0, tt/L, Fig. [6^,e. The 
corresponding solutions of the amplitude equation, Fig. 
|6lD-d,f-g, agree with the envelopes extracted from the di- 
rect solutions of the Helmholtz equation. Fig. [3l The 
eigenvalues of Eq. (IM)) also give the frequencies that 
correspond to the bandedge states that are also in excel- 
lent agreement (the observed deviation is less then 0.1%), 
Fig. [3l Knowlendge of the bandedge frequencies allows 
to determine all parameters of the tight-binding approx- 
imation for w(fc), Eq. ^. Therefore, the entire band 
structure in the spectral region of the flat bands can be 
obtained solely from the solutions of the amplitude-phase 
equation. 

Filled circles in panels a,e of Fig. [6] denote the spa- 
tial position where the particular 4){x) turns to rm: /2. 
At these special points dA{x)/dx = denoted by the 
vertical dashed lines in b-d and f-g panels. Thus, the 
overall phase accumulated by 0(a;) over one period, is an 
important parameter indicative of the spatial structure 
of the amplitude. At the bandedge frequencies of the 



i{N) 



j{N) 



bands A\^ ' and B^^ ' (see Section 3 for notations), the 
phase is a bound function \4){x) — 0(0) | < 7r/2. There- 
fore, a; = 0, i/2, L are the only positions where the corre- 
sponding amplitude function takes minimum/maximum 
values. Fig. [6)3, f, - A{X) has only one "hump" for A\ 
and B\ . Comparison of 0(x) computed for fc = (solid 
lines) and k — tt/L (dashed lines) edges of each photonic 
band shows (Fig. [6]) that the difference occurs in the 
spacial regions where the electromagnetic waves propa- 
gate via tunneling mechanism in language of CROWs of 
Section 3. In these regions A{x) is small, which explains 
the small spectral width of the corresponding photonic 
bands. 



As the eigen-frequencies of the higher order states 
A2 , B2 ... shift further away from the primary band- 
gap region, Fig. [TJd, 0(a;) becomes progressively steeper 
function, leading to the steady increase of the number 
of "humps" in A(x),[ni Such progression accelerates the 
spacial dependence of amplitude and eventual breakdown 
of the scale separation approximation used in the deriva- 
tion of Eqs. (|24l25p . Nevertheless, such loss of appli- 
cability occurs well outside of spectral region of interest, 
demonstrating the robustness of the developed approach. 



SUMMARY 

In this work we studied the optical properties of dual- 
periodic photonic superlattice with four different theoret- 
ical methods. Although each method has its limitations, 
the results obtained in framework of each model comple- 
ment each other. 

Numerical simulation with transfer matrices yields the 
direct result for the photonic band structure. This 
method, however, provides little physical insight into the 
nature of the photonic bands. In Section 3 we com- 
pared the spectrum of the infinite (periodic) crystal with 
the transmission spectrum of the finite system with the 
length equal to one period of the superstructure. We 
identified the individual transmission resonances with the 
photonic bands and found one-to-one correspondence. 
Furthermore, the spatial distribution of the fields at a 
resonance demonstrated that in L :^ a limit the enve- 
lope (amplitude) of the state changes slowly - on the 
scale of L. 

With a method commonly employed in condensed mat- 
ter physics we investigated the resonant interactions be- 
tween Bloch waves when the second, longer-scale, mod- 
ulation is introduced. We showed that flattening of the 
photonic bands is related to, but goes beyond band fold- 
ing. It arises due to the increasing coupling between 
Bloch waves with k- vectors at the boundaries of Brillouin 
zone. The subsequent increase of the band-gaps regions 
"squeezes" the the bands making them progressively flat- 
ter as A^ = L/a is increased. Although, as showed in the 
Section 4, this approach fails for very large N, it still pro- 
vided an insight into the origin of the anomalously small 
dispersion in spectra of PhSC. 

Diffraction gratings introduced in optical fibers are of- 
ten spatially modulated. Coupled-mode theory has been 
developed to reduce the problem to a study of the am- 
plitudes of the forward and backward propagating waves 
and to avoid the direct solution of Maxwell equations. 
Although the refractive index contrast induced in opti- 
cal fibers is orders of magnitude smaller then in the sys- 
tems we are concerned with in this work, CMT-based 
approach of Section 5 provided a clear physical picture. 
It showed that the electromagnetic states of our optical 
resonators can be thought of as the eigen-states of the 



photonic wells. This further reinforced the analogy with 
CROWs that we developed in Section 3. 

We proceeded by noting a formal similarity between 
the Helmholtz equation with the considered dielectric 
function and the equation describing parametric reso- 
nance in oscillation theory. Adopting amplitude-phase 
formalism accompanied by separation of scales (short a 
and long L = N x a) allowed us to derive to tractable 
set of equations on the envelope functions. This enabled 
the study physically meaningful mode profiles directly, 
without assuming smallness of modulations of refractive 
index. 

We find that dual-periodic photonic lattice is a promis- 
ing system for studying slow-light phenomena. It can be 
viewed as an array of evanescently-coupled optical cav- 
ities. The cavities are formed due to presence of the 
second, long-range, modulation which creates alternat- 
ing spatial regions of allowed and forbidden propagation. 
Thus, analogies with multiple-quantum- well and coupled 
resonant optical arrays (CROWs) are appropriate. Un- 
like other CROW implementations, the studies structures 
can be produced with holography-based approach, which 
enables fabrication of arrays with a large number of iden- 
tical resonators. The latter is crucial condition from an 
experimental point of view. It should ensure that the 
resonances of the individual cavities couple to form flat 
photonic bands. 
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